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Abstract The quasi-periodic oscillations (QPOs) in black hole (BH) systems of different 
scales are interpreted based on the magnetic reconnection of the large-scale magnetic fields 
generated by the toroidal electric currents flowing in the inner region of accretion disk, where 
the current density is assumed to be proportional to the mass density of the accreting plasma. 
The magnetic connection (MC) is taken into account in resolving the dynamic equations 
of the accretion disk, in which the MC between the inner and outer disk regions, the MC 
between the plunging region and the disk, and the MC between the BH horizon and the disk 
are involved. It turns out that the single QPO frequency of several BH systems of different 
scales can be fitted by invoking the magnetic reconnection due to the MC between the inner 
and outer regions of the disk, where the BH binaries XTE J1859H-226, XTE J 1650-500 and 
GRS 1915H-105 and the massive BHs in NGC 5408 X-1 and RE J1034H-396 are included. In 
addition, the X-ray spectra corresponding to the QPOs are fitted for these sources based on 
the typical disk-corona model. 

Key words: accretion, accretion disks — black hole physics — magnetic fields — stars; 
individual (XTE J1859H-226, XTE J 1650-500, GRS 1915H-105) — galaxies: individual (NGC 
5408, RE J1034H-396) 

1 INTRODUCTION 

As is well known. X-ray quasi-periodic oscillation (QPO) is a common phenomenon in the radiation from 
black-hole binaries (BHBs). The high-frequency QPOs (HFQPOs) have been observed in several BHBs, 
some of which show interesting 3:2 frequency pairs (GRS 1915H-105, GRO J1655-40, XTE J1550-564, 
H1743-322), and others display single QPO (XTE J1859H-226, XTE J1650-500). QPOs have also been 
observed in ultraluminous X-ray sources (ULXs), e.g., a 54 mHz QPO in M82 X-1 and a 20 mHz QPO 
in NGC 5408 X-1 were discovered respectively by Strohmayer & Mushotsky ( 120031 ) and Strohmayer et al. 
(|2007t with XMM-Newton. Strohmayer & Mushotsky ( |2009t discovered another strong 10 mHz QPO in 
NGC 5408 X-1, and found that the correlation of timing and spectral properties of this source is similar to 
those of Galactic BHBs. 

The first convincing QPO of active galactic nuclei (AGNs) was discovered by Gierlinski et al. (I2008I I 
in narrow line Seyfert 1 RE J1034H-396, which opened a window for the comparative timing studies of 
stellar-mass and supermassive black holes (BHs). The 1-hour X-ray QPO observed in RE J1034H-396 is 
analogous to the 67 Hz QPO in BHB GRS 1915H-105 (Middelton & Done l20T0] i. In earlier years, quasi- 
periodic signals like QPOs were discovered in some supermassive BHs. For example, the power density 
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spectra of two X-ray flares of Sgr A* observed in 2000 and 2002 show five distinct peaks at periods of 100, 
219, 700, 1 150, 2250 seconds (Aschenbach et al. |2004i) and a quasi-periodic flux modulation with a period 
of 22.2 minutes was discovered in the X-ray data of the Sgr A* flare in 2004 (Blanger et al. 12006b . 

A number of models have been proposed to interpret HFQPOs in BHBs. However, none of them 
can fully explain the characteristics of QPOs, especially the correlations of spectral and timing properties 
(Remillard & McClintock l2006l Maiti-a & Miner l20T0] l. As is well known, HFQPOs in BHBs are strongly 
correlated to steep power law (SPL) state. A successful model of QPOs should link the corresponding spec- 
tral state, allowing for a highly dynamical interplay between thermal and nonthermal processes with the 
mechanisms operating over a wide range of luminosity (Remillard & McClintock |2006l l. QPOs observed in 
massive BHs are probably related to the bright state which is similar to the SPL state of BHBs (Strohmayer 
& Mushotskv 120031 120091 Middelton et al. 120091 ). and they may have the same origin as the HFQPOs in 
BHBs (Gieriinski et al. |2008l Bian & Huang [20T0l ). 

Inverse Compton scattering is generally thought to be the promising radiation mechanism of SPL state 
(Zdziarski 120001 McClintock & Remillard 20061 hereafter MR06), and disk-corona model is successful in 
fitting the spectrum of SPL state (Gan et al. 120091 hereafter G09; Huang et al. l20T0l l. Zhao et al. ( 120091 ) 
(hereafter Z09) interpreted HFQPOs in BHBs as the magnetic reconnection of large-scale magnetic fields 
generated by toroidal electric currents flowing in the disk without considering the spectral state. In this 
paper, we improve the model of Z09 based on the disk-corona model with the magnetic connection (MC) 
effects. We fit both the QPO frequencies and the corresponding X-ray spectra of four BH systems of dif- 
ferent scales, in which two BHBs (XTE J1859H-226, XTE J1650-500), one ULX (NGC 5408 X-1) and one 
supermassive BH (RE J1034H-396) are included. 

This paper is organized as follows. A description of the model is given in Section 2. The QPO fre- 
quencies and X-ray spectra of the four sources are fitted based on the disk-corona model with MC effects 
in Section 3. And some issues related to our model are discussed in Section 4. Throughout this paper, the 
geometric units G = c = 1 are used. 

2 MODEL DESCRIPTION 

2.1 Origin of magnetic fields in BH systems 

Large-scale magnetic fields play very important roles in high energy astrophysical phenomena. The rela- 
tivistic jets from AGNs and BHBs are launched and collimated by invoking the open large-scale magnetic 
fields related to the BZ and BP processes (Blandford & Znaiek lT9771 Blandford & Payne [1982]), and the 
broad iron lines observed in BH systems could be interpreted by transferring energy and angular momentum 
from a spinning BH to its surrounding accretion disk by invoking the magnetic connection via the closed 
large-scale magnetic fields connecting the horizon with the inner disk (Wilms et al. 120011 Miller et al. 120021 
Li I2002al Wang et al. 120021 hereafter W02). However, a consensus on the origin of large-scale magnetic 
fields in BH systems has not been reached. 

Unlike neutron stars, magnetic fields cannot exist on the horizon of an isolated BH, and they should 
be maintained by the surrounding environment, such as an accretion disk. For BHBs, the magnetic fields 
probably come from the plasma of the companion. One of the most promising origins of large-scale mag- 
netic fields in BH systems is accretion disk around the BHs, and some authors calculated the magnetic field 
configurations by assuming the toroidal electric current in the disk (Li l2002bl Z09). However, the origin of 
the currents remains to be clarified. 

In this paper, we intend to interpret the origin of the electric current based on the assumption that the 
accreting plasma deviates somehow from electric neutrality as it flows through the Lagrange point into the 
Roche lobe of the BHBs, resulting in toroidal electric currents flowing in the accretion disk. In addition, 
we resolve the dynamical equations of the accretion disk by taking the MC effects into account and figure 
out the mass density and current density in the disk, and then determine the configuration of the large-scale 
magnetic fields by considering the interaction between the electric current and the disk with the iterative 
algorithm. It turns out that we can fit the association of QPOs with spectral states in BH systems of different 
scales based on our model. 
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Fig. 1 The configuration of large-scale magnetic 
fields generated by the continuous distribution 
of toroidal electric current flowing in the disk. 
The magnetic field lines are plotted in thin lines 
and the boundaries between three types of mag- 
netic connection (MCHD, MCPD, and MCDD) 
are shown by the thick lines. The red line repre- 
sents the magnetic field line connecting the radii 
Tin and ro. The figure is plotted with meH = 10, 
a, ^ 0.8, m 0.1, a = 0.3, r/ = 10^"^ and 
n = 5, where rii and a are the mass accretion 
rate in unit of the Eddington accretion rate and 
the viscosity parameter, respectively. 




Disk 



Black Hole 



Fig. 2 A schematic drawing of magnetic recon- 
nection in MCDD, where the red line represents 
the field line connecting radii Hn and rg . 



The deviation from electric neutrality is described by defining a parameter 77 as follows, 

77 = |ne - npl/np, (1) 

where rie and rip are the number densities of electrons and protons, respectively. The charge density can be 
expressed as 

Pe = verip, (2) 

where e = 4.8 x lO^^^e.s.u is the electron charge. 

The large-scale magnetic fields generated by the toroidal electric current may in turn affect the current. 
For example, the field lines may pipe hot electrons into the corona above the disk therefore change the 
charge density in the disk. This effect should be strongest in the inner disk because the magnetic field 
intensity is strongest near the inner edge of the disk and decreases rapidly outwards. For simplicity, we 
assume this effect decreases with the increasing disk radius as a power law and the surface density of the 
current in the disk can be expressed as 

J = Pe • 2/i • i^K • r-" = r,enp ■ 2h ■ Vt^r^-"^ = ^p^hftKr^-'', (3) 

/imp 

where h, and nip are the half height of the disk, the mass density and the proton mass respectively. The 
power law index n is a free parameter for fitting the QPO frequency, and /i = 0.615 is the weight-average 
molecular weight of the gas. The toroidal electric current is assumed to distribute from the inner edge of 
the disk to the outer boundary of the disk-corona system, the radius of which is set at rout — lOOM in 
calculations, since the radiation of the accretion disk mainly comes from the inner region. 
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Fig.3 The contour of 190 Hz QPO of XTE J1859+226 in 77 — n parameter space with ?t!bh = 12, 
a* = 0.96, m = 0.38, and a = 0.3. 



Following the work of Znajek (|1978l l. Linet (1 19791 1. Li (|2002bl l and Z09, we can calculate the toroidal 
component of the electric vector potential determined by the current given by Eq. (|3) with a given mass den- 
sity of the disk matter And the boundaries of inner-outer regions related to the magnetic field configuration 
can be determined. As shown by Fig.l, there are three types of magnetic field configuration generated — 
MC of the BH with the disk (MCHD), MC of the plunging region with the disk (MCPD), and MC of the 
inner and outer disk regions (MCDD). In the MCDD region, the magnetic field lines are frozen in the disk 
at the inner and outer footpoints. The field lines will twist themselves since the inner and outer footpoints 
have different angular velocities. In Fig. 1, the inner and outer footpoints of the field line in red locate at Hn 
(the inner edge of the disk) and ro, respectively. The value of the angular velocity difference between the 
footpoints of this field line is maximal among all the lines in MCDD region, so this line will twist itself first 
and trigger the magnetic reconnection, as shown in Fig. 2, which releases magnetic energy periodically to 
generate flares with the frequency of the difference between Keplerian frequencies of the inner and outer 
footpoints. We interpret this frequency as the QPO frequency, and it reads: 



vqpo 



^K(nn) - ^K(n)) 
27r 



(4) 



where 



rigj[3.23 X 10^'^Hz and ttibh = M/Mq is the black hole mass in unit of the solar mass, 
a, = a/M is the dimensionless BH spin. r\n is initialized at the innermost stable circular orbit (ISCO). 

As shown by Eq. (O, parameters n and rj determine the current density and therefore magnetic con- 
figuration and QPO frequency once the surface density of disk matter is given. The value of 77 determines 
the current intensity and therefore magnetic field strength, while that of ?i determines the concentricity of 
the currents and magnetic fields. A larger n leads to a smaller MCDD region, a smaller distance btween rin 
and ro, and hence a lower QPO frequency, while a smaller n leads to the opposite results. As an example, 
a contour of 190 Hz QPO of XTE J1859H-226 is plotted in 77 — ?i parameter space as shown in Fig. 3. It is 
found that a maximum jy ~ 10^^^ is required to avoid a too strong magnetic field for a stable solution with 
the MC effect. 



2.2 The dynamical equations of disk and corona with MC effects 

The dynamical equations of the accretion flow are modified by considering the MC effects on the transfer 
of energy and angular momentum (Li |2002al G09), 



dr 



dr 



(5) 
(6) 
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Fig. 4 An equivalent circuit for MCDD, MCPD and MCHD given from left to right. 



where Q and g are respectively the viscously dissipated energy per unit disk surface and interior viscous 
torque of the disk, and Q consists of two parts, 

Q = Qd + Qcon (7) 

where Qd is radiated from disk as black body, 

Qd = fTTjf, (8) 

and Qcor is transported into the corona to heat it by magnetic reconnection of tangled small-scale magnetic 
fields, and it reads (Liu et al. |2002| l 

Q - - (9) 

where Bu, Va and p are the intensity of the small-scale magnetic fields, the Alfven speed, and mass density, 
respectively. 

The quantity H^c in Eqs. (5) and (6) is the flux of angular momentum and it reads 

1 dTMc 1 /d*\' n,-n,+i 



Hmc ^ — = — — —, (10) 

Anr dr 2'Krdr \2tt J dZi 

where dTMc and d^* are the torque exerted on to an infinitesimal annulus of the disk and magnetic flux 
threading the infinitesimal annulus, respectively. 

We can calculate a variety of magnetic transfer of energy in BH accretion disk by using an equivalent 
circuit in an analogous way to Macdonald & Thorne (|1982| l and W02 as shown in Fig.4. Three types of MC 
(MCDD, MCPD and MCHD) indicated in Fig. 1 are given in Fig.4, in which a series of loops correspond to 
the adjacent magnetic surfaces arising from the rotation of the closed field lines. 

The quantities dZu, dZp, dZa and dZc are the resistances corresponding to the two adjacent magnetic 
surfaces in the disk, the plunging region, the BH horizon and the corona, respectively. The quantities deo 
and dee are respectively the electromotive forces generated by the rotation of the disk and the BH horizon 
as given by W02, and dep is the electromotive force due to the rotation of the plunging region. 

The quantities fli, and dZi in Eqs. (6) and (10) are explained for the three different disk regions 
as follows. 

(1) MCDD region: Hi and ili+i are respectively the angular velocities of the inner footpoint i and 
outer footpoint i + 1 of the closed field line on the disk, and the resistance dZj is defined by 

dz,=dz,,, = :^, (11) 

where R^or is the average areal resistivity of the disk and corona, which is related to Ru, the surface resis- 
tivity of the BH horizon as follows, 

i?H = 47r = 377ohm, ^ ^ 
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where ryR is a parameter to adjust the value of i?coi- For simplicity, we set jjr = 0.1 in calculations. 

(2) MCHD region: $7, = and ili+i = are the angular velocities of the BH horizon and the disk 
respectively. The resistance 

dZ, = dZn + dZeor, (13) 

where 

dZn = 2pnde/w. (14) 

The parameters in Eq. (14) are given by 

w = (Eh/ph) sin 9, p'^ = r'^ + cos^ 

Sh = 2A/rH, re = M{\ + q),q= ^l-al (15) 

(3) MCPD region: fli = ftp and fi^+i = ilo are respectively the angular velocities of the plunging 
region and the disk, and flp is given as follows (Shapiro & Teukolsky[l983 Wang |2000l l. 

{r-2M)L^ + 2aME^ 
P {r^ + a^r + 2Ma^)E^ -2aMLr ^ ' 

The resistance is the same as that in MCDD expressed by Eq. (11). 

We adopt the same assumption as given by G09 about the corona: the optical depth of the corona 
Tcoi = 1 and the height of the corona / = lOrms. 



3 FITTING THE QPO FREQUENCIES AND X-RAY SPECTRA 

In our model, the toroidal electric current interacts with the dynamics of the accretion disk. To solve the 
dynamic Eqs. (5) and (6), we must know the configuration of the large-scale magnetic fields generated by 
the toroidal electric current, whose distribution is related to the surface density of the disk matter by Eq. (|3), 
and pni is in turn figured out by solving the dynamical equations. For simplicity, we assume that initially 
there is no electric current in the disk, and solve the Eqs. (5) and (6) with Hy^c = to get the global solution 
of the disk-corona system. The emerged spectrum of the disk-corona system is then simulated with Monte 
Carlo method based on the code developed by G09. The free parameters of the disk-corona model, e.g., 
BH spin a, and mass accretion rate m, can be determined by fitting the observed spectrum. So the surface 
density of the disk matter is obtained. Then we consider the interaction between the electric current and 
disk-corona with the iterative algorithm which consists of the following steps: 

(i) Assuming a value of the power law index n in Eq. (|3]l, e.g. n = 5. 

(ii) Calculating the surface density of the toroidal electric current in the disk and the configuration of 
the large-scale magnetic fields. 

(iii) Resolving the disk-corona system by taking the MC effects into account. 

(iv) Repeating steps (ii) and (iii) until the surface density of disk matter remains stable. 

(v) Calculating the frequency vqpq by Eq. dUi. 

(vi) Repeating steps (i)-(v) until the QPO frequency is in accordance with the observations. 

Steps (ii)-(iii) should be repeated several times before the surface density of disk matter becomes stable. 
Finally, we simulate the emerged spectrum again until it is unchanged. 

To avoid the negative radiation flux from the inner disk due to the transfer of energy and angular 
momentum to the outer disk by the MCDD as argued by Gan et al. ( |2007| l. we adjust the radius of the inner 
boundary of the disk to deviate outwards from ISCO. Although the deviation is less than 10% of the radius 
of ISCO, it results in ~^30% decrease of QPO frequency. 



3.1 Fitting the HFQPOs in XTE J1859-H226 and XTE J1650-500 with SPL state 

Observations show that HFQPOs in BHBs are associated with the SPL state which is characterized by high 
luminosity (therefore high accretion rate), strong power-law component and the steep power-law index 
(r > 2.4). We fit the single HFQPOs and the corresponding X-ray spectra of two BHBs XTE J1859H-226 
and XTE J 1650-500. The comparisons between the simulated spectra and the observed ones are shown 
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Fig. 5 The simulated spectra of XTE J1859+226 with different parameters: (a) tobh ~ 12, 
a* = 0.96, m = 0.38 and a = 0.3; (b) mBH 12, a* = 0.998, m = 0.27 and a = 0.3; 
(c) rriBH = 7.6, a* = 0.998, m ~ 0.31, and a = 0.3. The total spectrum and its thermal 
and Comptonized components are plotted in solid, dotted and dashed lines, respectively. The 
observation data are taken from MR06. The source distance is set at 11 kpc (Zurita et al. 120021 
MR06) and the inclination i — 70° is assumed. 



for XTE J1859+226 and XTE J1650— 500 in Figs. 5 and 6, respectively. The total spectra consist of the 
thermal component emitted from the disk and the power-law component generated by the inverse Compton 
scattering of the soft photons by the relativistic electrons in the corona. 

The fitting parameters are listed in Table 1, in which m and a are respectively the mass accretion rate in 
terms of the Eddington accretion rate (1.4 x lO^^meH g s^^) and the viscosity parameter The BH mass of 
XTE J1859H-226 was estimated in the range 7.6-12AfQ (MR06). We fit the spectrum of XTE J1859H-226 
observed in October 16— 18 of 1999 (MR06) with the BH mass of the upper limit 12il/0 and of the lower 
limit 1.6Mq. Since the BH mass of XTE J1650— 500 has not been well estimated, we take a larger mass 
7.3Mq and a smaller one 4-AIq, which are estimated by Orosz et al. (I2004l i to fit the spectrum. Very large 
BH spin is needed to fit the spectra for the lower limit to the BH mass, and the lower limit to BH spin 
corresponding to each BH mass is given in Table 1. The spectrum of XTE J1650— 500 is fitted with the 
lower and upper limits to spin for each BH mass, while the spectrum of XTE J1859H-226 can be fitted only 
for an extreme Kerr BH with the lower limit to the BH mass 7.6Mq. The value of the viscosity parameter 
a is selected properly in the range 0.1—0.3 to fit the spectrum of each source. The difference of the value 
arises probably from the difference of the strength of the small-scale magnetic field in the disk if the viscous 
process is dominated by the tangled small-scale magnetic field. The values of the fitting parameters n and 
T] for QPO frequencies are listed in the last two column in Table 1 . We fix 77 at a certain value for each 
scale of BH mass. BH systems with smaller mass have stronger magnetic field allowing of larger value of 
rj. A larger value of n is needed with a smaller BH mass, implying that the toroidal electric currents are 
concentrated very close to ISCO. 

3.2 QPOs in NGC 5408 X-1 and RE J1034-h396 

One of attractive features of our model is that it is applicable to fit QPOs with corresponding X-ray spectra 
observed in the BH systems of different scales. Strohmayeret al. ( 120071 ) discovered a strong 20 mHz QPO 
in the ULX NGC 5408 X-1, and the X-ray timing and spectral properties are analogous to those exhibited 
by Galactic stellar-mass BHs in the 'very high' or SPL state. The fitting parameters for the QPO and X-ray 
spectrum of NGC 5408 X- 1 are listed also in Table 1 . The simulated spectra for different parameters and the 
comparisons to the observed energy spectrum are shown in Fig. 7. However, at present, the BH masses in 
the ULXs are constrained to a wide range from 10'^ Mq to IO'^Mq (Miller et al. 120031 Cropper et al. 120041 
Roberts et al. |2005l Wu & Gu |2008] l. We take a very large BH mass lO^Mg to fit the data of NGC 5408 X-1 
since both the X-ray spectrum and QPO frequency are fitted better for larger BH mass. The QPO frequency 
and spectrum are both taken from Strohmayer et al. (12007) with the same Hydrogen column density, which 
are fitted in our model with the lower and upper limits to the BH spin as shown in Table 1 . 
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Fig. 6 The simulated spectra of XTE J1650— 500 with different parameters: (a) tobh = 7.3, 
a* = 0.87, TO = 0.10 and a = 0.2; (b) tobh = 7.3, a* = 0.998, m = 0.03 and a = 0.2; (c) 
Tobh = 4, a* = 0.91, to = 0.17, and a = 0.2; (d) tobh = 4, = 0.998, to = 0.06, and 
a ~ 0.2. The observation data are taken from Miniutti et al. (|2004t . The source distance is set at 
2.6 kpc (Homan et al. 120061 1 and the inclination i ~ 70° is assumed. 



The first convi ncing QPO of AGNs was reported by GierUnski et al. ( |2008l ) in RE J 1034+396. 
Middelton & Done ( 1201 Oi l suggests that the QPO discovered in RE J 1034+396 has an analogy to the 67 Hz 
QPO seen in the BHB GRS 1915+105 due to their similar 'hot disk dominated' energy spectra. Unlike other 
HFQPOs, the 67 Hz QPO in GRS 1915+105 is an exceptional case, which appears in thermal-dominant 
(TD) state (MR06). In this subsection, we fit both the 67 Hz QPO in GRS 1915+105 and the 0.00027Hz 
QPO in RE J 1034+396 as a c ompari son. The spectrum of GRS 1915+105 showing the 67 Hz QPO is 
taken from Middleton & Do ne (l20T0l l. And the spectrum of RE J1034+396 showing the 0.00027Hz QPO is 
taken from Middleton et al. (I2009I I with the same minimum galactic absorption fixed at 1.31 x 10^°cm^^. 
The 'hot disk dominated' spectra of GRS 1915+105 and RE J1034+396 are both fitted with almost max- 
imum BH spins as shown in Table 1, and their comparisons with the observed spectra are shown in Figs 
8 and 9, respectively. The disparities between the simulated spectra and the observed ones in the energy 
bands 10-30 keV for GRS 1915+105 and 0.5-1 keV for RE J1034+396 indicate probably that a second 
Comptonization process is needed to generate the TD spectrum as shown in Figs 10 and 1 1 of Middelton & 
Done (12010b . where a low temperature, optically thick thermal Comptonization is added to fit the spectra. 
This Comptonization may be generated from the transition layer between the disk and the corona. 

Inspecting Table 1 and Figs 5—9, we find that the QPOs and the corresponding X-ray spectra of BH 
systems of different scales can be fitted with the magnetic reconnection of the large-scale magnetic fields 
based on the disk-corona model. QPOs in massive BHs have similar features with those of BHBs, e.g. very 
centralized distribution of the electric currents or large value of n, and its inverse proportion to BH mass. 
They are either related to the SPL state, like XTE J1859+226 and XTE J1650-500, or related to the TD 
state, like GRS 1915+105. These similar characteristics hint probably that QPOs in BH systems of different 
scales may have the same origin and are associated with the same spectral state. 
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Fig. 7 The simulated spectra of NGC 5408 X-1 with different parameters: (a) tobh = 10^, 
a* = 0.95, m = 0.012 and a = 0.1; (b) meH = 10^ a* = 0.998, m = 0.005 and a = .1. Th e 
plot style is the same as Figure 5. The observation data are taken from Strohmayer et al. (12007V 
The source distance is set at 4.8 Mpc (Karachentsev et al. |2002| l and the inclination i = 75° is 
assumed. The total Hydrogen column density is set as nn = 13 x 10^° cm"^ (Strohmayer et al. 
I2OO7I 1. 



4 DISCUSSION 

In this paper, a toy model for the QPOs in BH systems of different scales is proposed based on the magnetic 
reconnection of large-scale magnetic fields generated by the toroidal electric currents in the disk. The dy- 
namical equations of accretion disk are resolved based on the interaction between the electric currents with 
the disk-corona system by using the iterative algorithm. The 190 Hz and 250 Hz single HFQPOs in BHBs 
XTE J1859+226 and XTE J1650— 500 associated with SPL states are well fitted based on the disk-corona 
model with electric currents flowing in the inner disk. And the similar QPOs observed in ULX NGC 5408 
X-1 and Seyfert 1 AGN RE J1034H-396 and the corresponding X-ray spectra are also fitted. The spectrum 
of NGC 5408 X-1 is fitted with strong power-low component and steep power-law index suggesting that the 
QPO is similar to the HFQPOs in BHBs XTE J1859+226 and XTE J1650-500 and are probably associated 
with the same spectral state — SPL state. While the QPO in RE J1034+396 is analogous to the 67 Hz QPO 
in GRS 1915+105 which is associated with TD state. 
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Fig. 8 The simulated spectra of GRS 1915+105 with different parameters: (a) tobh = 18, a* = 
0.998, m = 0.25 and a = 0.14; (b) tobh = 10, a* = 0.998, m = 0.4 and a = 0.16. T he plo t 
style is the same as Figure 5. The observation data are taken from Middleton & Done (12010b . 
The source distance is set at 1 1 kpc (McClintock et al. I2006I I and the inclination i = 66° (Fender 
etal. [T999l lis used. 
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Fig. 9 The simulated spectra of RE J1034+396 with different parameters: (a) mBH = 7 x 10®, 
a* = 0.99, m 0.15 and a = 0.3; (b) ttibh = 7 x 10^ a* = 0.998, m = 0.14 and a = 0.3; (c) 
TTiBH = 2 X 10®, a* = 0.998, m = 0.13, and a = 0.3. The plot style is the same as Figure 5. The 
observation data are taken from Middleton et al. (12009) . The source distance is set at 125.9 Mpc 
(0=0.042, 7?o=100 km s~^ Mpc~^) and the inclination i = 40° is assumed. We only consider 
the minimum galactic absorption fixed at 1.31 x 10^° cm~^ (Middleton et al. 120091 ). 



Similarities of the QPOs in BH systems of different scales enable us to estimate some physical quan- 
tities of the massive BHs which have not been well constrained at present. As suggested by MR06, the 
frequencies of HFQPOs in BHBs showing the 3:2 frequency pairs are inverse proportional to the BH mass. 
Abramowicz et al. (|2007| l used the 1 /M scaling to expect QPO frequencies for BHs of different scales and 
neutron stars. Similarly, we can use this relation to estimate the BH mass in ULX NGC 5408 X-1. Including 
the three BHBs showing the single HFQPO, we have the relationship, jtibh = 1460t'Qpo, between the BH 
mass and the QPO frequency as shown in Fig. 10. The BH mass of NGC 5408 X-1 is then estimated about 

7.3 X lO"*Af0 with the 0.02 Hz QPO as shown by the black dot in Fig. 10. 

The quasi-periodic signals similar to QPOs were also discovered in the X-ray flux of Sgr A*. If these 
signals are indeed QPOs and are triggered by the magnetic reconnection described in this paper, then the 
lower limit of the BH spin can be constrained because the outer footpoints of the magnetic field lines cannot 
extend to the infinite distance. We estimate the lower limit of the BH spin of Sgr A* as 0.448 by fitting the 
22.2 minutes signals discovered in the X-ray flare on 2004 August 31 (Blanger et al. 12006b using the mass 

4.4 X 10® Me (Genzel et al. l20T0l ). which is very close to the value a* k, 0.44 ± 0.08 estimated by Kato et 
al. ( 12010b using the QPO method in the context of disk-seismology. 
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It is noticed that the upper and lower kHz QPOs of accreting X-ray binaries are interpreted as the 
Alfven wave oscillations with different accreted material mass densities at a preferred radius near the star 
surface, and this model successfully explains the empirical relation between the upper and lower kHz QPO 
frequencies and the linear relation between the high and low QPO frequencies of BHs, neutron stars and 
white dwarfs (Zhang [20041 Zhang et al. |2007l l. The idea of interpreting QPOs with magnetic reconnection 
proposed in this paper may also apply to neutron stars and other astrophysical objects. For example, QPOs 
can be interpreted as the reconnection of magnetic field on the surface of a neutron star, then the model 
evolves to the beat frequency model (e.g.. Miller et al. 119981 ), providing a physical mechanism for QPO 
production. 

Magnetic reconnection is being increasingly recognized as an important process in high energy objects, 
such as stellar X-ray flares, accretion disk corona, and magnetar flares. In this paper, we apply it to interpret 
the QPOs in BH systems of different scales based on a possible origin of the large-scale magnetic field 
in BH accretion disk. The magnetic fields generated in the inner disk by the electric currents are ~ 10^ 
and ~ lO'' Gauss for BHBs and AGNs respectively. Although the strength of these fields is much smaller 
than that of the tangled small-scale magnetic fields in the disk, the magnetic reconnection of the large-scale 
magnetic fields should have some effects on the heating of the corona, and these should be considered in 
the future work. Furthermore, the influence of the magnetic reconnection on X-ray spectra is another open 
question to be solved. 
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